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Abstract 

The concept of quasiresonance was introduced in connection with inelastic collisions be- 
tween one atom and a vibro-rotationally excited diatomic molecule. In its original form, the 
collisions induce quasiresonant transfer of energy between the internal degrees of freedom 
of the diatom: there is a surprisingly accurate low order rational value for the ratio of the 
changes in the vibrational and rotational classical actions, provided the vibrational and rota- 
tional frequencies of the diatom are approximately related by low order rational values, and 
the collision was longer that the rotational period of the molecule. In this paper we show 
that quasiresonance can be extended to many other processes and systems, and that it may 
be understood in terms of the adiabatic invariance theory and the method of averaging. 

1 Introduction 

The term quasiresonant vibration-rotation energy transfer (QRVR) was first introduced p][2] in 
the context of low velocity inelastic collision processes between one atom and a vibro-rotationally 
excited diatomic molecule. These are transitions which dominate all other transitions over a 
significant range of initial conditions and occur when there is a highly specific and efficient 
inelastic transfer of energy between the internal vibrational and rotational degrees of freedom 
of the diatom. 

The first experimental evidence of QRVR was reported in a level-resolved study of vi- 
brotationally inelastic collisions between Li^ and noble gas atoms. The experimental results 
show remarkably sharp, strong peaks in the vibro-rotational inelastic rate constant at very spe- 
cific high initial rotational states with values that substantially exceed the purely rotationally 
inelastic rate constant. A extremely good correlation between the change in the rotational (J) 
and vibrational (v) quantum numbers of the diatomic at the most important peak of the final 
rotational states was also observed, with a very strong propensity rule Aj = — 4Av. 

Since then QRVR transfer has been observed in a wide variety of collision species and inter- 
action potentials at different energy ranges [31 13 • It has been suggested as a possible relaxation 
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mechanism of rotationally excited light molecules, such as H2 molecules formed in laboratory 
traps or interstellar clouds, and good candidate to provide a rotational inversion mechanism in 
rotationally pumped lasers at ultracold temperatures [Q. Also, a similar quasiresonance energy 
transfer mechanism has been shown to control the process of predissociation for the most weakly 
bound state of van der Waals complexes jS] . 

In general, the qualitative features that characterize the quasiresonant vibration-rotation 
energy transfer are a large vibrationally inelastic cross section with very narrow final rotational 
state distributions, a strong correlation 

Aj = - Av, (1) 

with p/q a rational fraction involving small integers. The correlation becomes stronger at higher 
initial rotational states and lower collision velocities. Finally, experimental results obtained for 
different collision species also indicate a remarkable insensitivity to both the initial vibrational 
state and the nature of the interaction potential ^ EH [1] . 

These striking properties of the quasiresonant vibro-rotational energy transfer effect have 
stimulated a number of theoretical studies [21 El HI El E3 ED3 ITT] , both classical trajectories 
analysis and quantum mechanical computations, of inelastic cross sections and rate coefficients 
for collisions between one atom and a rotationally excited diatom at ordinary and ultracold 
temperatures. The good qualitative agreement between classical and quantum results suggested 
that the classical dynamics of the system is playing a significant role in the mechanism underlying 
the QRVR energy transfer. 

Classical trajectory studies [21 El El have shed light on the dynamics involved in this process. 
It has been observed [2] that collisional processes that occur at high molecular rotation states and 
low collision velocities are composed of a sequence of small collisions, or collisionettes, in which 
the interaction potential presents remarkable peaks. These collisionettes take place whenever 
the rapidly rotating diatom and the external atom come closer to each other, with the diatom 
in its outer vibrational turning point and the molecular axis nearly collinear with the atom. 
In between collisionettes the interaction potential becomes insignificant and the rotational and 
vibrational actions remain nearly constant. 

The ratio between the change in the two molecular actions is approximately (and often 
spectacularly close to) a low order rational ratio of the vibrational and rotational classical 
molecular frequencies. This suggests that the vibrational and rotational frequencies of the diatom 
are in approximate low-order resonance when that QRVR transfer occurs, and consequently a 
commensurate relation of the type 

q rij ujj 

holds, with n v and rij small integers and u v and Uj the classical vibrational and rotational 
molecular frequencies. 

The observed frequency-locking relation (J2J) suggested that the interaction potential is act- 
ing as a nonlinear external coupling between the internal vibrational and rotational molecular 
oscillators. This conjecture was later supported by Hoving et al. [H] who explained QRVR 
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transfer in terms of adiabatic invariance and presented numerical evidence which showed that 
this phenomenon is associated with the appearance of very large isolated nonlinear resonances 
in the phase of the system. Numerical results also indicated that the general rule (pQ) followed by 
the molecular rotational and vibrational quantum numbers in the dominant QRVR transitions, 
which according to © can be expressed as 

n v Av + rijAv = (3) 

is connected to the conservation of the action 

I = n v v + rijj (4) 

Quasiresonant energy transfer processes have also been observed in atom-diatom collisions 
that occur at ultralow energies [Z|. However, in this regime there are new aspects that introduce 
qualitative differences between the classical and quantum behaviors. The strong correlation 
between Aj and Av persists, but the magnitude can be much smaller than one quantum, making 
the quantum transitions classically forbidden. In the quantum limit, significant threshold effects 
appear in the quenching rate coefficient in the zero temperature limit causing some of the most 
relevant quasiresonant channels to remain closed due to energy conservation. Classically, the 
collisionettes described in processes at higher energies are replaced by a strong modulation of 
the interaction potential at the characteristic frequency of the quasiresonant transition. For 
certain initial conditions the incoming atom can get temporarily trapped around the diatom, 
bouncing several times before being rejected by the repulsive component of the interaction 
potential. Although the temporal evolution of the atom along these trapping trajectories is a 
extremely sensitive function of the initial conditions (chaotic scattering region), remarkably the 
quasiresonant correlation (^Q) is still satisfied. 

In this work we will focus on the classical analysis of quasiresonance. We will show that this 
is a common effect, not restricted to atom-diatom inelastic collisions, which arises from tran- 
sient "turning on" of interactions between quasi-resonant, i.e. not necessarily exactly resonant, 
degrees of freedom. The same quasiresonance effect will be shown to arise from an explicitly 
time dependent transient "turning on" and one that occurs autonomously. 

Although the work presented bears some relation with adiabatic switching the con- 
nection has not been explored. Adiabatic switching is a method for semiclassical quantization 
which had some success in starting with a multidimensional separable problem of known quan- 
tized actions and switching on interactions adiabatically, hoping that the actions do not change, 
and thus allowing the energy with those actions to be read off of the trajectories which have 
suffered the switching. 

The paper is organized as follows. We begin in section |2] with the theoretical analysis of 
classical system under a transient interaction. In section we illustrate the main features of the 
quasiresonance effect numerically, using a two-oscillator resonant Hamiltonian, a new example 
of quasiresonant system. Section |1] further analyzes the mechanisms behind quasiresonance, 
especially from the point of view of phase space. Here we also present a more detailed discussion 
of the quasiresonance effect between the vibrational and rotational internal degrees of freedom 
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of a diatom molecule in slow inelastic collisions with an external atom. Finally, in section El we 
summarize the main conclusions of the work. 

2 Adiabatic invariants under a slow transient interaction 

In this section we show that quasiresonance is a common effect which arises from the cou- 
pling of approximately resonant internal degrees of freedom of a classical system, perturbed by 
some transient interaction with additional degrees of freedom. We start by considering a two- 
dimensional integrable system Hq, perturbed by the transient interaction with one structureless 
particle. The Hamiltonian of the total system may be expressed as 

H = H (3) + H k (P) + eV(J,Q,P,Q&) (5) 

where J = ( J\, J2) and <!> = ($1, ^2) are the action-angle variables that describe the unperturbed 
system Hq, Q and P the coordinate and the conjugate momentum (not necessarily action-angle 
variables) of the incoming particle and its kinetic energy, e is a small parameter which 
characterizes the magnitude of the interaction term. We will assume that the term satisfies the 
condition 

ton V(3,Q,Q&) =0 (6) 

^— >oo 

with t, = i;(J, <I>, P, Q) an interaction parameter that controls the amplitude of the interaction 
between the system Hq and the particle. In a typical collisional process this parameter would 
correspond to the distance between centers of mass of the colliding species, which do not notice 
the interaction potential unless they are close enough to each other. 

We are interested in the analysis of a collisional event in which the transient interaction 
of the system Hq with the particle induces a change in its internal state from the initial state 
(Ji,3>i) = (Jii, J2i,*ii,*2i) at t = Ufci -» 00) to the final state (J f ,<S> f ) = (Ji/, J 2 f, $i/$2/) 
at t = tf (<^f — » 00), see figure Q. 

Let us assume that the two independent frequencies of the unperturbed Hamiltonian Hq 
satisfy an approximate resonant condition of the form 

Moj 2 ~ Nlox ~ (7) 

where 

f) H 

^ (J) = ib? ( * = 1 ' 2) (8) 

The secularity in the unperturbed Hamiltonian due to this resonance condition can be removed 
applying standard secular perturbation theory |13j . Namely, we can perform a canonical trans- 
formation from the variables (3,<f>) to a new set of action-angle variables (I, (|)), which define a 
frame that rotates with the resonant frequency. Taking the generating function for the canonical 
transformation as 

F(I,P,*,Q) = (M$ 2 - iV*i)Ii + $2^2 + P ■ Q (9) 
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t = t, t = tf 

Figure 1: A sketch of the initial (a) and final (b) states in the transient interaction process between the 
system Hq a structureless incoming particle. 

the equations for the transformation between the two set of variables of the Hamiltonian Hq are 
given by 

dF 

fa = — - = M$ 2 - N$i (10) 

Oil 

dF 

^ = wr ^ (11) 
dF 

J > = a¥T - Nh <12) 

dF 

J 2 = l^ = M/ l+ / 2 (13) 

In terms of the new variables the transformed Hamiltonian may be expressed as 

H' = flo(I) + A*(P) + eV(I,^P,Q;& (14) 

If the approximate resonance condition (JJJ) is satisfied, the evolution of the new angle variables 
in the rotating frame is given by 

/ 8V dV\ 
fa = M$ 2 - AT$! ~ e (M— - iV— J (15) 

and 

dV 

fa = $ 2 = w 2 (I) + e— (16) 

Hence, provided that the frequency uji is far enough from zero, the oscillation of the new angle 
variable fa near the resonance will be fast compared to the variation of the angle variable fa. 
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Under these conditions, the dynamics of the system in the proximity of the nonlinear resonance 
can be described by the averaged Hamiltonian 

H = flo(I) + H k (P) + eV(I, l5 P, Q; (17) 

where 

0i, P, Q; = / * <|), P, Q; £)d0 2 (18) 

Since H does not depend on the angle variable (j>2 its conjugated action variable 1% is an adibatic 
invariant of the motion in the proximity of the resonance. But, according to 1)12(1 and (|13|). this 
adiabatic invariance of the action conjugated to the fast angle variable <f>2 implies that 

M 

h = J2 + j-jrJi = const (19) 

or 

NAJ 2 + MA Ji = iV ( J 2f - hi) + M ( J lf - J u ) = (20) 

which gives the rational ratio in the action changes that is observed in the quasiresonance effect. 
Thus, a quasiresonance is implied by the existence of an adiabatic invariant characterizing the 
quasiperiodic motion of the system, in the proximity of a nonlinear resonance between its internal 
degrees of freedom. 

From this perspective quasiresonance is merely the survival of an adiabatic action, corre- 
sponding to a fast angle. However, subtleties regarding the boundaries of adiabaticity as probed 
by the strength and duration of the interaction which couples the unperturbed actions leave a 
rich subject to be explored. 



2.1 Example: Atom-vibrotor collision 



A nice illustration of quasiresonance effect can be found in the dynamics of an atom, fixed on 
a spring and is rotating on a plane, which is slightly perturbed by the collision with a slow 
incoming particle. Before the interaction occurs the dynamics of the unperturbed rotating atom 
is described by the integrable Hamiltonian 

* = *.<*,*) = £ + 5^3 +PM (21) 

where m is the atom mass, p r and pe are the momenta conjugated to the radial and angular 
variables r and 8, U(r) the internal elastic interaction potential, and 



If, 1 



„2 



Jr = ^zf Prdr = — j>\\ 2m [H - U(r)] - ^dr (22) 



2-7T / 2n J V T 

and 



Jd = 2vr / Pedd = Pd ^ 
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a) 




Figure 2: In a) the trajectory (a fixed "gear" shape) described by an atom fixed on a spring when there 
is a integer ratio of six between the vibrational and rotational frequencies. This is the form of the fast 
angle-averaged potential in the adiabatic analysis. In b) is a sketch of the evolution of the vib-rotor after 
it is perturbed by a slow colliding particle. 



the vibrational and rotational action variables. 

Initially the atom is vibrating in and out as it rotates on a plane, tracing out a "gear" 
shape which is actually fixed in space if the radial and angular frequencies satisfy a resonance 
condition (0), see panel a) in figure @. More generally, the "gear" is slowly rotating clockwise 
or counterclockwise initially, with angular frequency <fi±, i.e. the slow angle. The gear-shaped 
potential is just the potential of the averaged Hamiltonian given in Eq. 1171 

The interaction with the incoming particle induces a coupling between the internal degrees 
of freedom of the spinning atom and, therefore, a change in its vibrational and rotational fre- 
quencies. In the proximity of a resonance, when the condition Muj r — Nojg ~ is satisfied, 
such variation produces a slow rotation of the gear, which starts moving as an effective rigid 
rotor, see panel b) in figure ©. The angular momentum associated with this rotating gear is 
given by the new action variable 1$, while the action conjugated to the rapid angle variable in 
the rotating coordinate system, in this system the radial or vibrational action, is the modified 
adiabatic invariant that satisfies the condition (|19[) . Namely, 



I r = J r + ^ Je = ^ j J 2m [H - U(r)] - ^dr + ^-p 9 = const (24) 

This example illustrates the physical nature of the reduced dimensionality system which results 
when the average over the fast variable is performed. As long as the quasiresonance holds, the 
reduced dimensionality object acts consistently as if it did not possess the degrees of freedom 
removed by averaging. 
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2.2 Nonautonomous systems 



The physical picture of the "gear" as the averaged Hamiltonian strongly suggests that an external 
perturbation pushing non-reactively (nonautonomously) could play the role of the incident atom 
in the example above, causing the gear to change precession in the same way as did collision 
with the atom. Certainly the canonical transformation to slow and fast angle variables is still 
relevant, as that involves only the "system". Thus we consider an integrable system perturbed 
by an explicitly time dependent interaction, i.e. the interaction parameter is i; = t,(t)- In terms 
of the action-angle variables (J,*!*) of the unperturbed system, the total Hamiltonian may be 
expressed as 

H = H (3) + eV(J,Q;%) (25) 

where e is a small parameter that characterizes the perturbation strength, which we will assume 
that satisfies the transient condition 

lim V(I,<S>;Z,(t)) = (26) 

t— >±oo 

In order to remove the singularity associated with the resonance (J2J) between the unperturbed 
frequencies Q, we again seek a canonical transformation from the variables (J, O) to a new vari- 
ables (I, (j)) in the rotating coordinate system. We choose for such transformation the generating 
function 

F(I, = (Af$ 2 - iV$i)ii + $ 2 / 2 + /(I, <£, (27) 
where / is a multiply periodic function of the angles that can be expanded as the Fourier series 

/(I, a nin2 (I, Qeiim*!-***') (28) 

711,712 

With such generating function, the equations of the transformation between the two sets of 
variables are 

6 X = — = M$2 - iV$i + (29) 

^ = I = ** + f < 30 > 

■ h =wr ~ Nh + si <31> 
J * = m 2 = Mh+h+ lk (32> 

In terms of the new variables, the transformed Hamiltonian can be written as 

H' = H(l-X) + ^ = H(m+^ (33) 
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where ^ = a%/dt. In the proximity of the resonance, the evolution of the new angle variables is 
given by 

• d 2 f ( dV dV\ ■ d 2 f 
* = M* 2 - ^ + ^ K c - N-) + ^ (34) 

and 

^ 2 = * 2= „ 2(I)+£ |K + ^ (35) 

Let us now assume that the transient interaction is a function with slow time variation and, 
therefore, we can consider ^ as a small parameter that characterizes the slow change in the 
interaction. Provided that the frequency 0J2 is sufficiently far from zero, (|34|) and 1|35|) indicate 
that the angle <f>2 oscillates much faster than <f>i near the resonance. Under these conditions, the 
slow motion of the system in this region can be described by the averaged Hamiltonian 

Thus, the application of the averaging over the fast angle leads to its conjugated action as 
the first term of the series for the modified adibatic invariant of the system near a nonlinear 
resonance. As in the autonomous system, such a modified invariant reflects the change that the 
transient perturbation induces in the topology of the phase space trajectories in the proximity of 
a nonlinear resonance zone. To zero order in the perturbation strength and slowness parameters, 
the propensity rule associated with a quasiresonance arises also in a non-autonomous transient 
interaction from the linear relation between the modified invariant and the action variables of 
the unperturbed system. 



3 Numerical studies and further analysis of quasiresonance 

We have just seen that the quasiresonance can arise from a explicitly time dependent or an 
autonomous transient interaction. A well defined model system suitable for further study is 
composed of two anharmonic Morse type oscillators that are transiently coupled by a non- 
autonomous external perturbation. In terms of the action-angle variables of the two uncoupled 
oscillators, we express this model resonant Hamiltonian as 

H = H (J) + V(<t>,t) (37) 

with the Hamiltonian for the unperturbed system 

H Q (3) = a\J\ + C12J2 - anJf - CL22J2 (38) 

where en, an > (i = 1,2), and the non-autonomous transient perturbation as a sum of resonant 
coupling terms 

m m 

V(&,t) = v g(t) £ sin(n 2 $ 2 - n x $i) (39) 

ni=l ri2=l 
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with ri2, n\ and m integers and vq the coupling strength. We will assume that the coupling 
interaction is turned on and turned off by means of a time-dependent a Gaussian function 



<?(£)= exp -(t-t p ) 2 /2a 2 



(40) 



characterized by a parameter a. 

Figure © displays the variation of the action J2 with respect to its initial value for trajec- 
tories with initial angle variables chosen at random. We observe that the action remains nearly 




10 20 30 40 50 1 2 3 4 5 6 



J 2 i % 

Figure 3: On the left, the variation of the action J2 (top) and the ratio of the changes in the two 
actions, AJ1/AJ2 (bottom), with respect to the initial value of the action J2. Both pictures include the 
data of 20 classical trajectories with the same initial actions and angle variables chosen at random. On 
the right, the J2 — 62 phase plane corresponding to the trajectories included in the panels on the left. 
The time-dependent interaction potential was frozen at time value t = t p . The total energy is E = 50, 
the Gaussian parameter a = 80, the coupling strength v Q = 0.05 and m = 3. The different parameters 
in the unperturbed Hamiltonian Hq are = 1 and an = 0.001 (i = 1, 2). All quantities are expressed in 
arbitrary units. 

constant except in a well defined region of initial action values where significant changes occur. 
This region corresponds to values of the initial action where the quasiresonance effect arises, 
with a low order rational value for the ratio of the changes in the actions of the unperturbed 
system. Here the correlation between the variations of the actions of the two oscillators in the 
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quasiresonance defines a large plateau region characterized by the propensity rule 

AJi + AJ 2 = (41) 

Figure (jHJ) also shows a view of the J2 — #2 phase plane when the coupling strength reaches 
its maximum amplitude at t = t p . From the analysis of the previous section it is clear that 
quasiresonance must be connected with the presence of resonance zones in the phase of the 
system [S], and here we see them explicitly. Although this surface of section was generated 
by examining the motion with the coupling artificially fixed at the maximum value, later we 
shall examine the temporal evolution in phase space under the actual transient interaction. As 
did Hoving and Parson [5] for atom-diatom collisions, we find that the region of phase space 
corresponding to the values of the initial actions where a quasiresonance occurs is dominated by 
a strong isolated nonlinear resonance at fixed coupling. Here, we observe a large resonance zone 
composed of three islands which are centered at the initial action value which gives the position 
of the center of the plateau region defined by the propensity rule (|4*Tj) . 

3.1 Resonance parallelogram 

Trajectories can get caught inside the limits of a growing island and thereby subject to action 
transport. Thus, the parallelogram shape structure observed in figure © for the change in 
the actions of trajectories included in a quasiresonance domain is an image of a nearly uniform 
spread of the final action values within the limits of former resonance islands, see figure @. 

3.2 How close is the action ratio to low order rationals? 

The hallmark of quasiresonance is the ratio of action changes adhering to a low order rational, 
reflective of a near (but not necessarily exact) resonance in the system phase space. How close 
to the rational value are the action ratios in fact? Linear plots are of no use, since the accuracy 
is typically far better than the linewidth of the plot. As figure © indicates, the significant 
variations in the actions for the initial conditions included in a plateau are correlated with 
relatively large changes in the energy of the system. Naively, variations of similar order of 
magnitude could be expected to occur in the ratio of action changes. However, even when the 
resonant condition is only approximate, the quasiresonant condition Mu)2 — Nuj\ ~ may be 
satisfied to much higher accuracy. The energy is also not conserved (in the non-autonomous 
systems) or it is exchanged between the perturber and the system (in the autonomous systems), 
but the action change ratio remains rational to a high degree of accuracy, often one part in 
10 6 for example in the present case. As seen in figure ©, the deviation of the ratio of action 
changes in a plateau from a low order rational value can be several orders of magnitude below 
the change in the energy or the individual actions. 

At higher energy the phase space which is accesible to the system is larger and new regions of 
initial action values that satisfy different quasiresonance conditions may appear. This is shown 
in figure ©, where five regions that have significant changes in the actions can be seen. Each 



11 



AJ 




! \ a.i=-.j.!.i_ 



J_ J c & J + J. 

Figure 4: On the left, a sketch of a resonance island in J— 8 phase space at constant interaction strength. 
On the right, the shaded parallelogram represents the region that confines all the possible values of the 
change in the actions, A J = Jj — Jj, for the trajectories with initial actions Ji within the limits of the 
resonance island on the left. The thick lines are the diagonals of the parallelogram, which are given by 
A J = and A J = -2 J, + 2J C . 

one of these regions defines a plateau in the ratio of action changes with a different rational 
value, r = N/M = 2/1, 3/2, 1/1, 2/3 and 1/2, in the propensity rule Q2UJI that characterizes 
the quasiresonance. Besides these main plateaus, a staircase structure composed of a series of 
small plateaus separated by regions of transition where the ratio of actions changes dramatically 
(quasiresonance is violated) is seen. As figure © also shows, these secondary plateaus are 
associated with quasiresonances characterized by large integer values M and N in the propensity 
rule (|2Uj) . Although in this figure the changes in the actions at these higher order quasiresonance 
domains are very small, and some of them are hardly visible, we will see below that they become 
more important as the coupling between the internal degrees of freedom of the system increases, 
and that they are involved in the eventual demise of the quasiresonance effect. 

As to be expected jB], the system with multiple quasiresonance plateaus now presents a 
more complex phase space structure, see figure (JZJ), Distinct resonance zones produce their own 
quasiresonance domains. Immersed in the phase space regions between the main resonance 
zones it is possible to distinguish some long chains of tiny islands associated with the secondary 
plateaus that appear in figure 
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Figure 5: From the top to the bottom, the ratio of action changes, the absolute value of the difference 
of the absolute value of this ratio from the rational number r = 1, the absolute value of the change in the 
action variable J2 and the absolute value of the change in the energy of the system vs the initial action 
Jn for one of the trajectories included in figure ©. 

3.3 Interaction strength and quasiresonance 

In the perturbative regime, or weak coupling limit, when the energy of the unperturbed inte- 
grable system is much higher that the strength of the transient interaction, each main quasires- 
onance domain remains isolated. The corresponding plateaus in the ratio of action changes are 
well defined and are punctuated by regions where the ratio of actions changes varies wildly. It 
should be noted that in such non-quasiresonant regions an approximate adiabatic invariance of 
both actions of the unperturbed system holds, but the ratio of whatever small changes do occur 
is not locked at rational ratios. 

As the strength of the coupling between the internal degrees of freedom of the perturbed 
system increases, the main quasiresonance domains grow, and higher order nonlinear resonances 
are activated in the regions between them. In the phase space, the generation of new quasires- 
onance domains is reflected by the destruction of tori where the actions of the unperturbed 
system are adiabatic invariants and the emergence of chains small islands of regular motion 
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Figure 6: The change in the action J 2 and the ratio of action changes vs the initial action J 2% for 
trajectories with the same parameters as in figure (jSJ), but energy E = 200. On the bottom, the histogram 
with the distribution of the values of the ratio of actions changes presented in the central panel. The 
dashed lines indicate the locations of some low order rational values r = n/m (n and m small integers), 
which coincide with the positions of the higher values in the histogram. 
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Figure 7: The J2 — #2 phase plane corresponding to the trajectories included in figure ©■ The time- 
dependent transient interaction was frozen at time value t = t p . 

associated with the modified adiabatic invariants in the proximity of higher order resonances. 
The combined processes of generation of new higher order quasiresonances and the expansion 
of the existing ones make the edges of neighboring quasiresonance domains to get closer and 
finally come into contact. Afterwards, the overlap between different quasiresonance domains de- 
teriorates them progressively and eventually is responsible for their disappearance. Hence, the 
Chirikov criterion |14j on the overlapping of neighboring resonances for the onset of chaos can 
be applied to estimate the destruction of a quasiresonance domain. Once such overlap occurs, 
all the invariant tori between the island of the resonance zones are destroyed and regions of 
intermingled chaotic motion surrounding the islands of regular motion become visible. In the 
strong coupling limit, when all the resonance zones have came into contact, nearly all the phase 
space appears filled up by these stochastic regions and only some isolated islands of the former 
main resonance zones may survive. Under these conditions all regions of regular quasiperiodic 
motion shrink and therefore all adiabatic invariants cease to exist. This picture of the demise 
of quasiresonance is approximate, however, since the timescale for the transient interaction can 
attenuate the effects of chaos; see Section @J 

In our model resonant Hamiltonian, the effects of a stronger coupling between the two 
oscillators on the quasiresonance can be analyzed taking larger amplitudes of the interaction 
strength vq. As figures (JHJ) and @ show, the increase in vq allows larger changes in the actions 
Ji (i = 1,2) around the resonant values. This causes an expansion of the resonances zones in 
phase space and makes more visible the long chains of small islands associated with higher order 
resonances. If the different resonance zones are well separated, an increase in the size of the 
corresponding quasiresonance domains might be expected. However, instead of that, figure 



15 




Figure 8: The ratio of action changes and the change in the action J 2 vs the initial action J 2 i for two 
different strengths v of the transient interaction. Panels on the left correspond to vq = 0.1 and panels 
on the right to vo = 0.4. In both cases the energy is E = 75, the Gaussian parameter in H40|) is a = 100 
and up to in — 3 coupling terms were included in the interaction term. The different parameters in 
the unperturbed Hamiltonian Hq are = 1 and an = 0.001 (i = 1,2). All quantities are expressed in 
arbitrary units. 

shows is a deterioration of the rational values for the ratio of action changes. The regions of 
stochastic behavior that emerge as a consequence of such overlap can be clearly observed in 
figure Q. 

An alternative way to change the coupling between the two oscillators s to include more 
coupling terms in the transient interaction, while keeping constant its norm in the angle variables. 
In this way we avoid the rapid expansion of the main resonance zones and introduce reduce the 
destruction of the main quasiresonance domains. In figures (Jlfljl and Ijlljl . for example, we set 
the norm as 

Nl = ^n = ^ (42) 

m 

where 




by taking a normalized interaction term given by 

m m 

V(0,t) = -g(t) V V sin(n 2 $ 2 - m$i) (44) 
m , , 

ni=l ri2=l 

As figure (|10|) shows, the increase in the number of coupling terms modifies the distribution 
of the final action values inside the main quasiresonance domains and induces higher order 
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Figure 9: The J2 — O2 phase space for the trajectories included in figure © at t = t p , when the transient 
interaction reachs its maximum amplitude. The panel on the left corresponds to vq = 0.1 and the one on 
the right to vq — 0.4. 

resonances, which give new secondary quasiresonance domains. Nonetheless, the emergence of 
higher order quasiresonances and the stronger overlap between neighboring domains eventually 
cause their destruction. The phase space of the system presents clear manifestations of these 
effects, see figure (|ll(l . with an increase in the number of islands inside the main resonance zones 
and the proliferation of chains of small islands associated with high-order resonances. Also thick 
stochastic layers arising from the overlap of neighboring resonance zones become clearly visible. 

4 Phase space mechanism of Quasiresonance 
4.1 Energy dependence of the resonance condition. 

The analysis of the phase space of the system has already shown the correspondence between the 
location of the center of a plateau in the ratio of action changes in a quasiresonance and the center 
of a nonlinear resonance zone: the center of the plateau corresponding to the propensity rule 
(|20l) coincides with the location of the rational unperturbed torus whose independent frequencies 
satisfy a resonant condition of the form ((7J). From now on we will denote as J^ M (i = 1, 2) the 
values of the actions of the unperturbed system that give the location of this resonant torus at 
energy E. That is, the solutions of the energy equation 




(45) 



and the resonant condition Q. 
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Figure 10: The ratio of action changes and the change in the action J2 vs the initial action Jn for the 
normalized transient interaction Q44|l . Panels on the left were obtained with m = 2 and panels on the 
right with m = 5. In both cases the energy is E = 75, the Gaussian parameter in l|4(J|l is o = 100 and 
the coupling strength vq — 0.4. The parameters in the unperturbed Hamiltonian are the same as in the 
previous figures. All quantities are expressed in arbitrary units. 




Figure 11: The J2 — 82 phase space for the trajectories included in figure l|10|) at time t = t p . The panel 
on the left corresponds to m = 2 and the one on the right to m = 5. 
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Figure 12: Location in the E — J2 plane of the lowest-order N : M resonant tori of unperturbed the 
Hamiltoninan (|38f) . The values taken for the parameters of the model are a; = 1 and an = 0.001 (i — 1,2). 
All quantities are expressed in arbitrary units. 

In the model resonant Hamiltonian introduced in the previous section, for example, the 
unperturbed Hamiltonian Ho defines an integrable system in which the actions of the two un- 
coupled oscillators, Ji (i = 1,2), are conserved and the angle variables vary linearly in time, 
$i = Ui(3)t + <I>i0j with Ui the unperturbed frequencies 

Ui(J) = -Qj- = a i- 2a H J i (* = !> 2 ) ( 46 ) 

The four-dimensional phase space of the unperturbed system is completely stratified into two- 
dimensional tori defined by angle coordinates (3>i, $2) and radii (J\, J2), with all periodic orbits 
restricted to lie on those invariant tori whose independent frequencies satisfy a commensurate 
relation. Thus, according to (|38|) and (J7J), the actions on the rational unperturbed torus N : M 
are related by: 

Figure p2|) shows the location in the E — J2 plane of the lowest-order unperturbed resonant tori 
in the range in which the unperturbed frequencies are positive, < Jj < J|" ax = (i = 1, 2) 



and < E < E max = 




As soon as the transient interaction is applied the unperturbed tori of the system begin 
to be deformed, but as we know from the Kolmogorov-Arnol'd-Moser (KAM) theorem, not 
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all of them in the same way. For sufficiently small coupling strength and Jacobians d(o/d3 
different from zero, this theorem guarantees that those tori bearing conditionally periodic motion 
with incommesurate frequencies continue to exist, being only slightly distorted. Rational tori 
whose independent frequencies satisfy a commensurate relation (JJJ) are grossly deformed or even 
completely destroyed. 



4.2 Span of a quasiresonance domain 

We analyze the phase space structures associated with the frozen, or constant transient interac- 
tion strength limit of the Hamiltonian. Although the interaction strength is always changing in 
quasiresonance, fixed strength phase space plots are still useful, in particular for estimating the 
span of a quasiresonance domain. 

We use our model resonant Hamiltonian (|37|). It is convenient for our purposes to introduce 
dimensionless variables. Considering the scale of time as 

r = * (48) 
a 

and introducing dimensionless action variables 

ji =4 ( i=i > 2 ) ( 49 ) 

ah 

with E the initial energy of the unperturbed system, we write the scaled Hamiltonian 

h=- (50) 

as 

m m 

h = siji + s 2 j2 - snjf - S22J2 + zg(j) Y sin(n 2 $2 - (51) 

ni=l 712=1 

where Si = aai, su = Ea 2 au (i = 1,2), e = vq/E and 

g{r) =exp[-(r-r p ) 2 /2] , (52) 

with t p = tp/a. We will assume below that the total energy of the unperturbed system is much 
higher than the maximum coupling strength, e « 1. In this limit the analysis of the dynamics 
of the system using first order classical perturbation theory will provide a good approach to 
describe the coupling process between the two oscillators. 

Let us now assume that the two frequencies of the unperturbed oscillators satisfy the resonant 
condition ((7J) . The removal of the resonant variables applying the canonical transformation given 
by the generating function 

F = (M$ 2 - iV$i)Ji + <S> 2 h (53) 
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leads to the Hamiltonian 



h = (s 2 M - Sl N)h + s 2 I 2 - (s n N A + s 22 M 2 )lf - 2s 22 Ml x l 2 - s 22 I^ + (54) 



mm r , n /T \ 



III III S J 

+eg{r) Yl sin { f n2 

ni=l n2=l 



H1 N 02 + iV^ 1 



The averaging of this transformed Hamiltonian over the rapid oscillation of the new angle variable 
(j>2 near the resonance zone gives the secular Hamiltonian 

m 

h = (s 2 M - Sl N)h + s 2 / 2 - (sn^V 2 + s 22 M 2 )I 2 - 2s 22 MhI 2 - s 22 l\ + eg(r) sin(fc^i) (55) 

k=l 

Thus the dynamics of the system in the proximity of the resonance ./V : M can be described by 
the equations of motion 

^P- = (s 2 M - Sl N) - 2(s n N 2 + s 22 M 2 )h - 2s 22 MI 2 (56) 

(IT 

dl m 

— = -eg{r) ^ k cos(fc<£i) (57) 
fc=i 

Our previous numerical results show that the length of the plateau regions in the ratio of action 
changes can be estimated from the maximum size of the nonlinear resonance zones. As was 
expected, see figure (|13j) . the phase space of the averaged Hamiltonian reproduces the main 
islands of the resonance zone, but not the chains of small islands associated with higher-order 
resonances and neither any stochastic region. 

The fixed points of the secular Hamiltonian occur at the action value 

T N:M = (g 2 M ~ siN) - 2s 22 MI 2 
1 2(s n N* + s 22 M*) 

which corresponds to the resonant actions J^ M (i = 1, 2) of the unperturbed Hamiltonian IjMHjl . 
In the weak coupling limit, the dynamics of the system in the proximity of the resonance actions 
can be described by the pendulum Hamiltonian 



h p end = -(snN 2 + s 22 M 2 )(Alf-- M ) 2 + eg{r) £ sin(^i) (59) 
with A/f :M =I X - 



k=l 

The 2m fixed points of h pen d occur at the resonance action I\ = , the m elliptic fixed 
points, Q < x < 4>\ 2 < ■■ < 4>\ m < 2ir, correspond to maximum values of the function 



= X>in(A:0i) ( 6 °) 



k=l 
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Figure 13: Poincare surfaces of section for the total Hamiltonian 151|) (on the left top) and for the 
averaged Hamiltonian l|55|l (on the right top) at time t — £ p , total energy E — 75, Gaussian parameter 
a = 100, coupling strength vq = 0.2 and m = 3. On the bottom the corresponding variations of the 
action J^. 



and m the hyperbolic fixed points, < 4>ii < 4>i2 < •• < 'Pim < 27r, to its minimum values. 
The curve of the separatrix that passes through the hyperbolic fixed point (t>ik(k = 1> ■■^ m ) is 
given by: 



\ 





/(0i) 




(sniV 2 + 


s 22 M 2 ) 



(k = l,...,m) 



(61) 



Figure (|14|) shows the function f{4>\) and the curves of the separatrices T^ :M (t, cj)\) for different 
values of m. The case m = 1, with only one hyperbolic fixed point, corresponds to a single 
pendulum Hamiltonian. For m > 1 there are an increasing number of hyperbolic fixed points 
which develop a structure composed of m encapsulated separatrices curves. The outermost 
separatrix is associated with the deepest well in the function f{<pi) and passes through to the 
hyperbolic fixed point (j)^ m . Whereas the innermost one is associated with the highest well and 
passes through the hyperbolic fixed point 4>ii- 

The maximum excursion of the action I\ along each separatrix curve T^ [M (t, <f>i) occurs at 



the elliptic fixed point 
be written as: 



Thus the maximum variation of this action at each instant r can 



AlS M (r) 



(k = l,...,m) 



(62) 



(sniV 2 + s 22 M 2 ) 

Returning to the initial action-angle variables to study the structure associated with the 
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Figure 14: The function f((f>i) i|f)U|) and the curves of the separatrices T^ :M (r, 0i) (fc = 1, ..,m) (fHT|> 
for the pendulum Hamiltonian h pen d (|59|l for several values of m. The values taken for the different 
parameters of the model were; N = M = 1, an = 022 = 0.001, vq = 0.1 and r = t p . 



resonance iV : M in the ji — (j'2 — $2) phase plane, the center of the resonance zone is located 
a t ii = Jxe 1 O2 = j§ E )• According to the transformation relations (|12j) - l|13j) and the equation 
(|BT))) there are a total number of 2mN (2mM) fixed points, < <£f 1 < : < <£f 2 < $1 2 < 

- < *i,mN < n,mN < 2^(0 < $1 >X < < $1 >2 < CD« 2 < ... < $| mM < $« mM < 27T). 

Hence the curves of the separatrices associated with the hyperbolic fixed points create a chain 
of N (M) identical units, each of them containing m encapsulated separatrices curves. Figure 
H15j) displays some of these chains. 

The maximum deviation of the action j\ (j'2) takes from its resonant value ji^ 1 (J2 £ ) 
along each separatrix curve is the same at all the elliptic fixed points $>f (Vp = 0, .., N — 1) 
[<3?2, P m+i (VP = 0, .., M— 1)]. Thus it follows from ([T2|) - (fT3|) and (|62|) that the maximum variation 
induced by the N : M resonance zone on the actions j\ and j'2 along the curve of the separatrix 
k at time r can be expressed as: 



and 



Ai 2 ^ M (r) 



= 2N\ 

max \ 




7(#,l) " 


(sniV 2 + s 22 M 2 ) 


= 2M\ 

max \ 






( SuN 2 + S22M 2) 



l,...,m) 



(A; = 1, ...,m) 



(63) 



(64) 



We define J iE f. (i = 1,2) as the two extreme values of the action Jj along the curve of the 
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Figure 15: Separatrices curves (161(1 in the J2 — $2 phase plane. The values taken for the different 
parameters of the model were: an = 022 = 0.001, vo = 0.1 and time r = t p . The dashed line gives the 
location of the resonance action J%e at energy E — 250 (arb. units). 



separatrix k when the resonance zone reachs its maximum size at time r = t p . That is, 

1 



T N:M± _ T N:M , 
J i,E,k — J i,E =■= 2 



a T N:M / 



1,2 ; k = l,.,m) (65) 



We can now compare our analytical prediction for the size and shape of a primary resonance 
zone with the length of the corresponding plateau in the ratio of action changes in the quasires- 
onance effect. Figures (|16|) and (|17|) display the quasiresonance domain associated with the 
nonlinear resonance iV : M = 1 : 1 for different number m of coupling terms. As was expected, 
the resonance zone presents a chain composed of mM islands, two in figure (|16|) and three in 
figure Q17|). associated with the mM elliptic and mM hyperbolic fixed points located at the reso- 
nance action J\\- The thick lines on the Poincare surface of sections represent the curves of the 
separatrices T^ :M (t p , cf)i) (k = 1, ...m) (|61|) when the transient interaction reachs its maximum 
amplitude at r = r p . There is good agreement between the Poincare surface of sections obtained 
from the numerical simulations resonance analysis for the main resonance zones. Thus, in the 
limit of weak coupling between the two oscillators (e << 1) the expressions (|63|) and (|64|) give 
a good estimate of the maximum spread of the actions across the phase space region dominated 



|17|) show, the maximum varia- 
AJ# M (tp)l (i = l,2) m, 



by the main resonance iV : M. Importantly, as figures (|16j) and 
tion of the action along the curve of the innermost separatrix, 
provides a good approximation for the length of the plateau in the ratio of action changes in 
the quasiresonance. 

The phase space between the innermost and outermost separatrices corresponds to deviations 
from the rational value of the action changes giving the quasiresonance domain. This is due 
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Figure 16: On the left is shown the Poincare surface section in the J2 — $2 phase plane. The thick 
solid lines are the curves of the separatrices l|61|) for N : M = 1 : 1, k — 1, 2 and t = t p . The dashed lines 
give the values of the actions J^^t for the different separatrices. In the central panel, the ratio of 

action changes, AJ1/AJ2, vs the initial action J2- In the panel on the right the variation of the action J2 
with respect to its initial value. In these last two pictures the data of 20 trajectories with initial angles 
chosen at random are represented for each value of the initial action J2- All data were obtained using 
the same parameter values: energy E — 100, coupling strength vq = 0.1, temporal Gaussian parameter 
a = 200, m = 2, ai = a-z = 1 and an = a22 = 0.001. All quantities are expressed in arbitrary units. 
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Figure 17: The same as figure <| 1 f i|i but for m = 3. 

to the development of higher order resonance zones. The effect of these secondary resonances 
increases when the coupling between the two oscillators becomes stronger and eventually makes 
all the plateau region to disappear. 

4.3 Effect of timescale of the transient interaction 

In the previous sections we have analyzed the phase space of the frozen Hamiltonian. Now, in 
order to understand the behavior of the actions when a quasiresonance occurs, we focus on the 
time evolution of the system as the external transient interaction is switched on and off. 

The quasiresonance effect belongs to the larger class of adiabatic phenomena. If the interac- 
tion turns on and off too fast, adiabaticity will be violated. The timescale is set by the nature of 
the interaction potential and how often it is "visited" by the trajectory of the system, which in 
turn depends on the initial actions; see Fig. 1181 The period of the visitation of the interaction 
potential must be much shorter than the time scale of the transient perturbation. The interac- 
tion potential may be localized to one or a few small regions in both angles, ($i,$2)j as it is 
in an atom-diatom collision for example. Or, it may be more extended, greatly affecting which 
resonance zones survive to cause quasiresonance plateaus. 

The longer the interaction timescale, the more high-order resonances and chaos can be 
manifested; see figure ()19|) . Thus, there is a balancing act set up between interaction time and 
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Figure 18: Trace of a slowly precessing, nearly 11:2 resonant vibrotor (compare Fig. (J2J), in which the 
trajectory visits each "tooth" of the gear only every other pass. If the resonance had been 11:3, it would 
be every 3rd pass, etc. This makes the quasiresonance condition easier to break, since the timescale for 
the transient interaction to now has to be slower by a factor of two compared to an 11:1 resonance. 

the strength and form of the interaction. By keeping the transient interaction long compared 
to the internal time scales associated with the low-order resonances, but short compared to 
the time required to complete the long periods involved in the high-order resonances, we can 
maximize the low order resonance plateaus. Similar issues are faced in the adiabatic switching 
method 

4.4 Time evolution of the actions in a quasiresonance domain 

Figure Q2UJI shows the time evolution of an ensemble of trajectories with initial angle variables 
chosen at random and four different values of the initial action J2 . The parameters used are the 
same as in figure lTol which shows the phase space and action ratio. In panel (a), where the initial 
action coincides with the resonance action J 2 ^ > we observe the spread of the action values along 
a series of branches corresponding to the oscillations of the trajectories around the elliptic fixed 
points when they remain trapped in the resonance zone. Once the transient coupling interaction 
starts decreasing and the resonance zone shrinks, the actions settle on different values around 
the resonance action and a diffusion in action of the original ensemble is finally observed. 

When the initial actions are between the center and the edges of the quasiresonance domain, 
J^E 1 — — J^El^i see P ane ls (b) and (c), there is a relatively smooth diffusion of the 
action until the strength of the transient interaction reachs its maximum amplitude. Afterwards 
the ensembles evolve into two horizontal distributions symmetrically arranged with respect to 
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Figure 19: The ratio of action changes and the change in the action J2 vs the initial action J-ii for 
transient interactions l|39Jl with different time dependence. The panels on the left correspond to a Gaussian 
function (|40Jl with a = 50, the panels on the center to a = 100 and panels on the right to a — 500. In 
all cases the energy is E = 75, up to m = 2 coupling terms have been included in the interaction term, 
which has a strength vq = 0.1. The different parameters in the unperturbed Hamiltonian are the same 
as in the previous figures. All quantities are expressed in arbitrary units. 
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Figure 20: The action J2 vs. time for an ensemble of 2000 trajectories with the same initial actions 
and initial angle variables $1 and $2 chosen at random. In (a) the initial action is J2i = 52.78 (J21 — 
J2100, center of the nonlinear resonance zone), in (b) is J hi = 48 (J^'loo < 'hi < ^100,1 > inside the 
quasiresonance domain), in (c) is J21 = 44.3 {J<n ~ J^ioo i~' one °f the edges of the quasiresonance 
domain) and in (d) the initial action is J2i — 42.5 {J2100/2 < < «^2ioo "L between the curves of the 
two separatrices) . The amplitude of the transient interaction is maximum at £ = 750. In the four panels 
the same parameters as in figure <|16|) were chosen. All quantities are expressed in arbitrary units. 
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Figure 21: The J 2 — $2 phase plane (black dots) for a ensemble of 30000 trajectories at different time 
values. All trajectories have initial angle variables chosen at random and the same initial actions as the 
panel (a) in figure H2U|) . The angle variable ranges from to 4tt. The gray dots represent the analytical 
curves of the separatrices (I61|) at the same time values. The transient interaction is zero in panel 1, it 
is increasing from panels 2 to 4, it reaches its maximum amplitude in panel 5 and it is decreasing from 
panels 6 to 9 where has practically disappeared. All the parameters of the model are the same as in figure 



the center of the resonance zone. The spread in final action decreases as the initial actions 
approach to the edges of the quasiresonance domain, J<n — J^ei^- ^ n this re gi° n the ensemble 
evolves into a bimodal distribution strongly peaked at the final action values J2/ — J2i and 
J2/ — —J2i + 2 J^g 1 , which correspond to the two linear branches that give the diagonals of the 
parallelogram depicted in figure ®. 

As panel (d) shows, the dispersion of the action is less important for initial values that range 
from the boundary of the innermost separatrix to the outermost separatrix, J 2 ki — — 
J 2 E m ■ Here the distribution of the action values defines a smooth structure that grows when 
the strength of the transient interaction is increasing and then shrinks as soon as such interaction 
starts decreasing. Thus, when the transient interaction disappears the action of nearly all the 
trajectories recover their initial actions. Trajectories with initial action even further away, 
outside the limit of the outermost separatrix curve, J 2 % < J^Em an< ^ > J^Em "> s ^ a y outside 
the phase space region swept by the resonance zone and their actions remain approximately 
constant throughout the process. 

Figures (|2 11123 \) give much insight into the behavior of the actions observed in figure (|20|) 
when the initial action is within a quasiresonance domain, J 2 ei~ — — ^2E x' ■ They show 
that crossing the separatrices associated with the nonlinear resonance zone is the underlying 
mechanism that makes possible the significant changes in the actions in the quasiresonance effect. 
The breakdown of adiabatic invariance due to the passages through the separatrices of nonlinear 
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resonance zones is a quite general phenomenon which has been observed to play an important 
role in the dynamics of systems that are perturbed by a slowly time-varying interaction; such 
as the dynamics of charged particles in electromagnetic fields in plasma physics, the motion of 
asteroids in celestial mechanics or the propagation of short radio waves in the ionosphere, see 
|17j and references therein. 

4.5 Phase space evolution 

Phase space pictures add insight to the mechanism of action changes and quasiresonance. Figure 
(|21jl displays the evolution in phase space of one ensemble of trajectories with initial actions 
equal to the resonance action, J2i = J^e 1 - Hence, the phase space points of the ensemble are 
already inside the oscillatory region of the resonance zone when it develops from the resonance 
action J^e- As the resonance zone grows in time the phase points cannot entirely escape from 
this region, but they can cross the innermost separatrix curve, getting trapped in the internal 
loops and encircling the elliptic fixed points, or they can remain confined in the region between 
the two separatrix curves. Figures (|22[) and (|23j) include ensembles of trajectories with initial 
actions between the center and the edges of the quasiresonance domain, J^ei" < < ^2Ei + - 
In both sequences the ensembles are initially unperturbed and therefore the action remains 
nearly constant. But as the resonance zone grows, it gradually deforms the invariant tori, 
J2 = Jn = const, that keep the trajectories confined in phase space. Once the initial ensemble 
reaches the outermost separatrix associated with the nonlinear resonance, and provided that the 
resonance zone is still growing in time, some phase space points pass through the separatrix and 
go into the oscillatory region of the nonlinear resonance zone, i.e. their dynamics changes from 
a precession to oscillation. As in the previous sequence, once the phase points are confined in 
the oscillatory region they can either remain in the transition regions between two separatrix 
curves or go inside the internal loops by passing through the innermost separatrix. 

It is only after the resonance zone starts shrinking that the trajectories trapped in the 
oscillatory region can recover the rotational dynamics by re-crossing the outermost separatrix. 
Afterwards there is no way they can pass through the resonance zone again. Since the trajectories 
can randomly emerge to any point inside the resonance zone, the final effect of crossing the 
separatrices is the diffusion in action of the original ensemble around the resonance action value. 
As the transient interaction goes away, the destroyed unperturbed tori, Jj = const (i = 1,2), are 
replaced by the modified invariant tori, I2 = const, that characterize the quasiresonance effect. 
These new invariant tori confine the trajectories across the phase space region previously swept 
by the separatrices of the nonlinear resonance zone. Since such region ranges from J 2 g j - to 

Jjg^ ', the length of the quasiresonance domain, important variations of the actions may occur 
after the transient interaction goes away. 

As we mentioned before, the nonadiabatic effects that arise in the proximity of the edges of 
the quasiresonance domains are associated with the progressive approach of the trajectories in 
this region to the boundary of the nonlinear resonance zone. As figure (j23j) shows, these are 
trajectories that intersect the separatrices very close to the hyperbolic fixed points, where the 
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Figure 24: Time evolution of two trajectories of the pendulum Hamiltonian l|59|l with different Gaussian 
parameters in the time-dependent function g(t); a = 1000 in panel (a) and a = 4000 in panel (b). Both 
trajectories had the same initial action, inside a quasiresonance domain. The dotted line gives the 
evolution of the separatrix. Note the closer approach of the trajectories to the separatrix in the case of 
the longer interaction time. 

motion of the system becomes extremely slow, and no adiabatic invariant remains. Trajectories 
with initial actions inside a well defined quasiresonance domain, see figure (|22jl. intersect the 
separatrices far enough from the hyperbolic points and do not generate significant nonadiabatic- 
ity These trajectories can go deeper into the oscillatory region of the resonance zone and fully 
perform the quasiperiodic motion around the elliptic fixed points that gives the quasiresonance 
effect. But, the more slowly the transient interaction proceeds, the more trajectories get closer 
to the hyperbolic fixed points when the crossings of the separatrices occur and, consequently, 
the more extended across the quasiresonance domain the nonadiabaticity effects manifest. 

The increasing attraction of the trajectories with initial action inside the quasiresonance 
domains towards the unstable fixed points is clearly observed in figure (|24|). which shows the 
evolution of two trajectories of the pendulum Hamiltonian (|59(1 with the same initial action but 
different slowness of the transient interaction. Initially the pendulum is rotating and the angle 
variable is a monotonically increasing function, but as the strength of the transient interaction 
increases the rotation slows down. The more slowly the transient interaction evolves the more 
the trajectory sticks to the hyperbolic fixed point. When the energy of the separatrix equals 
the energy of the pendulum, the trajectory crosses the separatrix and the transition to the 
oscillatory regime occurs. Afterwards the pendulum starts orbiting around the elliptic fixed 
points describing large loops that cover an extensive region in phase space. Later, when the 
strength of the transient interaction is decreasing, a new intersection of the separatrix occurs and 
the pendulum recovers its rotational motion. The more slowly the transient interaction proceeds 
the closer to the hyperbolic fixed point the trajectory emerges from the oscillatory region. Thus, 
in the limit of extremely slowly varying transient interaction, when the closest approach to 
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Figure 25: Evolution of two trajectories of the pendulum Hamiltonian (|59fl in a extremely slowly varying 
transient interaction, a = 10000. In (a) the outgoing trajectory from the oscillatory region follows the 
unstable manifold that gives the branch A J = —2 J,; + 2 J c in the bimodal distribution depicted in figure 
(J2J, and in (b) the trajectory evolves along unstable manifold that gives the branch A J = 0. Both 
trajectories have the same initial action, but different angle variable. The dotted line is the separatrix. 

the hyperbolic fixed points occurs, the outgoing trajectories from the oscillatory region of the 
resonance zone asymptote to the two unstable manifolds that emanate from such points, see 
figure (|25|). It is the evolution along such manifolds which gives the bimodal distribution in the 
final action values previously discussed. 

5 Conclusions 

We have shown that quasiresonance is a common effect, not restricted to vibro-rotationally 
inelastic collisions between an atom and a diatomic molecule, which may arise in any classical 
integrable system that is perturbed by a transient interaction with additional degrees of freedom 
which couples quasi (not necessarly exactly) resonant internal degrees of freedom. We have 
also shown that quasiresonance effects arise from either an autonomous or an explicitly time 
dependent transient interaction. 

The analysis of the variations of the actions and the energy of the system in a quasiresonance 
domain reveals that the ratio of action changes is paramount, adhering extremely closely to a 
rational ratio involving small integers, with the energy playing a secondary role. 

In the limit of weak and slowly varying transient interaction, the classical adiabatic invariance 
theory and the method of averaging have been used to establish the link between the presence of 
a quasiresonance domain and the existence of a topologically modified adiabatic invariant that 
characterizes the quasiperiodic motion of the system in the inner oscillatory region of a nonlinear 
resonance zone in phase space. Within a resonance zone, the rational value in the ratio of action 
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changes that gives a quasiresonance follows immediately from the linear integer combination 
that defines the modified adiabatic invariant, i.e. the action of the fast motion, in terms of the 
action variables of the unperturbed system. Outside the resonance zones there is practically no 
response from the system to the perturbative interaction and the adiabatic approximation is 
good for both unperturbed actions, but the ratio of the tiny changes is not remarkable. 

The generalized pendulum Hamiltonian obtained from the averaging over the fast motion 
in the proximity of a resonance zone has provided a good model to clarify the dynamics of the 
system when a quasiresonance occurs. The analysis of the phase space of such Hamiltonian 
has shown that the maximum size of the phase space region swept by the innermost separatrix 
of the nonlinear resonance in the course of time gives a good estimate for the length of the 
corresponding quasiresonance domain. This analysis has also shown that the change in the 
dynamics of the system due to the crossing of the separatrices as they move in phase space is 
the mechanism involved in the significant variations of the actions observed in a quasiresonance. 

A natural extension of the work presented in this paper is the generalization of the concept of 
quasiresonance to systems with more than two or three dimensions. To approach such systems, 
we can return to the physical picture of the averaged Hamiltonian, and the nature of the inter- 
action, where the averaging is over all fast variables which remain after identifying the slow ones 
(we do not assert that this is always possible). Whenever the transient interaction timescale 
(autonomous or not) is slow compared to the time between successive visits to the interaction 
region, the quasiresonance concept may survive the extension to many dimensions. There may 
remain more than one slow variable, in which case the nonlinear evolution will take place in the 
slow variable space, perhaps preserving none of their individual actions. Nonetheless all the fast 
actions would be conserved. 

Another extension of this work is its connection with possible models of Arnol'd diffusion 13 . 
The picture would be of random but quasiresonant travel along resonance surfaces in action 
space, interrupted by non-quasiresonant but extremely small action changes at the edges of 
quasiresonant plateaus. 
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